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1 Introduction 

Dynamical fermion simulations with Wilson-type fermions are now possible at small 
quark masses, large volumes [1-8] and small lattice spacings. With statistical precision 
reaching the (sub-)percent level, an important uncertainty which remains to be carefully 
controlled is due to the finite lattice spacing a. In the non-perturbatively 0(a)-improved 
theory, this issue has been investigated in some detail. On the one hand, significant a- 
effects have been found at lattice spacings of around a = 0.1 fm [9, 10]. On the other 
hand, both in the high precision computations of the scale dependence of coupling, 
quark masses and other composite operators [11-13] (see ref. [14] for a review) and in 
a recent scaling test [15], such effects were not visible. 

In this paper we analyse this apparent contradiction and find that the difference 
is of a simple kinematical origin, teaching us a more general lesson. In order to keep 
a-effects small, it is important to ensure that the general conditions necessary for the 
application of the Symanzik expansion of lattice observables in powers of a are well 
fulfilled. We discuss the general issue in Sect. [2] and then illustrate it in Sect. [3] on a 
particular example where precise non-perturbative results are available: a Schrodinger 
functional correlation function. The lessons we learn there allow us to perform an 
improved computation of the renormalization of the axial current in Sect. HI which we 
check further by considering modified renormalization conditions (Sect. 14. 2p before we 
conclude. 

2 The Symanzik expansion and kinematics 

Symanzik's effective theory [16-18] is the fundamental tool for analysing the approach 
of renormalized lattice observables, Oiat to their continuum limit Ocont- As an example 
consider the asymptotic expansion of a (space-momentum zero) correlation function 
(e.g. eq. ^M) 

C\i,t{xo, a) "^ Ccont(a:o) + aCi(xo) + o^ C2(xo) + . . . . 

The functions C depend in addition to the kinematical argument xq on the intrinsic 
scale A = Aqcd = 0(300 MeV) of the theory as well as the quark masses Mj. In the 
above formula the powers of the lattice spacing, are in principle modified by a further 
a-dependence in Ci{xo) — > Ci{xo,a). It is a consequence of the dependence of the 
couplings in the Symanzik effective Lagrangian on a as well as the dependence of the 
coefficients of the effective fields on a. However, since QCD is asymptotically free, these 
coefficients depend only logarithmically on the lattice spacing when it is small. Such 
logarithmic terms are of minor importance for our discussion; they are dropped here. 

2.1 Point-to-point correlators at short distances 

Let us assume a large volume and a "point-to-point" correlator, i.e. Ciat(a^O)O) = 
o'^ ^^(01(0)02(2;)) with local operators 01,2- We focus on the short distance regime, 



xq <C 1/A,xo <C 1/Mj and an 0(a)-improved theory where Ci = 0. Then xq provides 
the only diniensionful parameter. On purely dimensional ground the leading correction 
then becomes 

a^C2(xo) "^ const, x (a/xo)^ Ccont(2;o) • (2.1) 

It is enhanced at small xq. In other words what we have used here is that up to 
logarithmic corrections the short distance correlation functions have a unique power- 
law behaviour in xq. 

A particularly simple example is 

C[,^{xo,a) = -a'^{P-{0)P^{x)), with P"(x) = V^(x)75ir'^V'(x) , (2.2) 

X 

which one may, in principle, consider for renormalizing the pseudo-scalar density. One 
expects that the short distance behaviour C^^^ ~ const, x Xq is difficult to reproduce 
accurately on a lattice. In fact, one term in the Symanzik expansion originates from an 
O(a^) correction to the field 



P^ff (x) = P'^(x) + a^C2 dld,P-{x) + ... (2.3) 



which contributes 
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C^^{xo)a' ^^"^ 24c2^CZtM + -.. (2-4) 

Xg 

to eq. ()2.ip . An order of magnitude enhancement is due to the second derivative of the 
steep function Ccont- Even if C2 may be arising only at 1-loop of perturbation theor}o, 
the considered case suffices to illustrate our main point: cutoff effects may have a 
significant kinematical enhancement. Particular examples are correlation functions with 
strong short distance singularities. 

Of course, this is the reason why the connection between the perturbative short 
distance regime of QCD and the non-perturbative long distance one is carried out recur- 
sively in the strategy of our collaboration [14,19,20]. It is then possible to have a/xQ <^ 1 
and xqA <^ 1. But furthermore, by making use of Schrodinger functional boundary con- 
ditions [21], one may easily construct correlators with a weak time-dependence [22]. Let 
us discuss this relevant issue in some more detail. 

2.2 Schrodinger functional correlators 

We now assume a finite volume with Dirichlet boundary conditions at xq = and 
Xq = T as explained in refs. [18,21,23]. These allow for the definition of gauge invariant 
boundary fields, for instance 

^" = I^EEc(x)^5^"C(y) (2.5) 



^this means C2 ~ const. / log(aA) 



constructed from the boundary quark (and anti-quark) fields C (C) &t ^0 = 0- I^i space 
we use periodic boundary conditions with a phase 9, 

ij{x + kL) = e'^ ^p{x) , '^{x + kL) = e"*^ '^{x) . (2.6) 

Because of the zero momentum projection of the boundary fields, the correlator 

fl^{xo.a) = -^{O^P%x)) (2.7) 

has a smooth behaviour at short distances (and small L) of the form 

/P,,(xo) = 3e-2v^^-o/^ {1 + 0{f{L)] , (2.8) 

in QCD with 3 colors and for massless quarks. There is no kinematical enhancement 
of a-effects. Together with the smooth background field introduced for the definition of 
the running coupling [23,24], this explains the very small lattice spacing effects in the 
running coupling and running operators mentioned before. 

We may also discuss the behaviour at large distances. There a saturation by few 
intermediate states in the spectral decomposition (see ref. [25] for a discussion of w, E) 

fU^o) = const. ^u;„^e--«^"-(^--«)^™^ (2.9) 

n,m 

will give an accurate description of the correlation function. The energies E^ are the 
finite volume eigenvalues of the QCD Hamiltonian in the pion sector and E^ are the 
energies of states with vacuum quantum numbers. Their dependence on L and the 
lattice spacing is suppressed. If the kinematics, given hy 6,xq, L is such that effectively 
a few states with energies up to aE ^ 1 contribute, the xo-dependence may again be 
strong and cutoff effects may be enhanced. However, the enhancement will not be as 
large as in Sect. l2.ll since there is no singular dependence at small xq. In Sect. [3] we will 
see quantitatively how the behaviour changes as L and xq are increased starting from 
L,xo < 1/A. 

We remark that the foregoing discussion is of course not in contradiction to the 
perturbative behaviour. In the perturbative region, xq is small and many states con- 
tribute significantly. But asymptotic freedom implies that their coefficients Wnm are 
fine-tuned such as to produce the smooth behaviour of eq. (j2.8p . This is often called 
quark-hadron duality. 

2.3 Energies and matrix elements 

The most common application of Symanzik's effective theory is to energies and matrix 
elements, e.g. extracted from eq. (|2.9p at large xq. The expansion of such observables, 

Olat ~ Ocont + a Si + a^ S2 + . . . , (2.10) 

will be valid and accurate when the relevant energies, momenta and masses are small 
compared to a^^. For the 0(a)-improved theory, a first scaling test [15] showed that 
indeed the corrections to the continuum appear to be reasonably small when a < 0.1 fm. 



3 The Schrodinger functional correlator f^^ 




Figure 1: The correlator f^^ in the perturbative region, L < 0.5 fm. The continuous curve shows the 
tree-level behaviour, while data points are for (f{L/2) — 3.33 (or g^{L) ~ 5.5j with filled circles for 
L/a = 24, squares for L/a = 16 and triangle for L/a — 12. A dashed line shows the position of the data 
with f{L/2) = 1.50, L/a = 24. 

We now proceed to discuss numerical results for /p in different regions of L (or 
(f'{L)). Let us start from the short distance, weak coupling, region. In the non- 
perturbative computation of the scale-dependence of composite operators, such as P", 
our collaboration chose Schrodinger functional boundary conditions with 6 = 0.5, T = L 
and a vanishing background field [12]. From these simulations we have precise data for 
/j^^ in a range of the Schrodinger functional coupling 0.98 < g^iL) < 5.5 [11, 12]. In 
order to cancel the renormalization factor Zp, we normalize by /j^^. at the midpoint, 
xo = T/2 = L/2. 

At g^ = 0, eq. (j2.8p yields the slowly falling exponential plotted in Figure [1] as a 
continuous curve. The non-perturbative results at g'^{L) ~ 5.5 are next to this tree-level 
curve, even though in general perturbation theory is not very accurate at this coupling. 
For example the renormalization factor Zp is known to be a factor 2 away from one. 
We note in passing that the very close agreement with tree- level at xq/T > 1/2 is 
somewhat accidental, as e.g. at a weak coupling of (f'{L/2) = 1.50 the curve is some 
15 % away (dashed line). All in all, the smooth behaviour anticipated in our general 
discussion, is found to a remarkable degree in the whole range of L^0.5fm. Cutoff 
effects are very small. Even at the shortest distances, xq/o, the a^ effects are only a 
few percent at (f'{L) = 5.5. At weaker coupling, (f'{L/2) = 1.50, the a? effects are 




Figure 2: The correlator f^^^ for L ~ 0.8 fm, T « 1.8 fm, 9 — 0.5 is represented by squares. L/a = 
8, 12, 16 data increase at small xq. A dotted line indicates the location of 9 = 0,L/a = 8. Circles show 
the behaviour at L ^ 1.2 fm, T ~ 1.8 fm, 5 = 0. Their errors are of the order of the symbol size. 



invisible within our precision of k, 2% (the data are not plotted in order not to clutter 
the graph). As explained in the previous section, this behaviour of the Schrodinger 
functional correlators goes hand in hand with small cutoff effects in the step scaling 
functions which describe for example the running of coupling and quark mass [12]. 

In physical units we have thus far investigated the region L = T^O.Sfm. For L 
just somewhat larger, L = 0.8 fm, T = 9/4 L, and 9 = 0, we have previously observed 
significant a^ effects in our computation of Z\ [10]. They were prominent in the statisti- 
cally significant disconnected contributions to Za, which can be shown to vanish in the 
continuum limit. As a consequence Za determined with or without the disconnected 
diagrams differed by almost 15% at a ~ 0.1 fm. Although the disconnected diagrams 
vanish quickly as a is reduced, an unpleasantly large ambiguity remained at typical 
values of a. For details we refer to the quoted reference. 

We now turn to that same kinematics and subsequently to larger L. But first we 
point out that in this region the dependence of f^^{xo) / f^^^{T /2) on 6 is insignificant 
compared to the effects we will discuss. An explicit example is provided by comparing 
the dotted line of Figure [2] with the L/a = 8 data points. For our numerical demon- 
stration we will thus freely use data at available values of 6. 

We start with correlation functions /j^^. at the same parameters as in [10], except for 
that we remain with 6 = 0.5. The behaviour of f^^, plotted in FigureEJ differs drastically 
from Figured]- see the tree- level curve as a reference. The non-perturbative correlators 
drop steeply (note that we use a logarithmic scale) and follow the characteristics of 
the described intermediate regime, where neither the smooth perturbative behaviour is 



realized, nor a single intermediate state is dominating (the latter would be seen as a 
straight line in the figure). For the coarsest lattice the logarithmic slope (effective mass) 
at xo = T/2 is as high as 0.8 in lattice units. In this situation, we find indeed very 
significant lattice spacing effects. It is then not surprising that there are also large a^ 
effects in the form of the mentioned disconnected diagrams. 

For an improved determination of Za, we would like to choose a kinematical re- 
gion, where (1) the a^ effects are better suppressed and (2) the Schrodinger functional 
simulations can be done close to or in the chiral limit. The simulations can reach the 
chiral limit when the infrared cutoff on the spectrum of the Dirac operator is sufficiently 
large. Since in the Schrodinger functional the infrared cutoff is dominantly controlled 
by l/T, we do not want to increase T significantly compared to the previous T ~ 1.8 fm. 
On the other hand, the energies E'!^ are expected to be decreasing with L; in fact in 
the small L limit they scale as L^^. This leads us to consider a somewhat larger L, 
namely L ~ 1.2 fm. Also in such a situation, namely with T = 3/2 L, we have simula- 
tion results [26], generated for the determination of improvement coefficients h^, h^ — b-p 
as well as the renormalization of the bare mass rriq following ref. [27]. The circles in 
Figure [2] show ff^^ in this kinematical situatiorcl. The function has a much slower decay 
and indeed only moderate lattice artifacts. The energies of the states dominating the 
correlator appear to be significantly lower - at least the logarithmic slope never exceeds 
a value of 0.3 at the coarsest lattice spacing. 

This kinematics is a good starting point for a redetermination of the renormalization 
factor Za with reduced intrinsic a^ ambiguities. 

4 New determination of Z^ 

The renormalization condition of the axial current is obtained by considering an axial 
Ward identity in exactly the same way as it was done in ref. [10]. One starts from a 
Schrodinger functional correlation function of the axial current with two pseudo-scalar 
boundary operators, 

^"H = ^E^W^5ir"a;(x-y)C(y) (4.1) 

^'"M = ^EC'W75^r'^^(x-y)C'(y) (4.2) 

x,y 



and performs an axial rotation of the variables in a region around the axial current [28]. 
Using isospin symmetry and PCAC one obtains for vanishing quark mass the Ward 
identity 

Z\f\A{xQ,yo,^) = /i(w) + O(a^) , xo > yo , (4.3) 



^ At L/a = 12 the massless point is not reached, but the behaviour at the smallest mass, shown 
here, does not change significantly when the mass is increased. See Sect. |4] for details. 
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in terms of the correlation function 

/iA(xo,yo,^) = -^^6'^^V'^^(0"^H(Ai)S(x)(^i)^(y)0^(a;)), (4.4) 

/iH = -\{0'\u;)0'^{u:)). (4.5) 

The superscript I reminds us that the improved axial current 

{Ai)l = Al + acAd'^P\ ^^(a;)=^(x)7^75^rXx), (4.6) 

is to be inserted. We refer to section 2 in ref. [10] for a derivation, the exact lattice 
implementation as well as a generalization to finite mass, which we use. 

Compared to that work we only change the kinematics. First of all the computation 
is now performed on a lattice of size L ~ 1.2 fm instead of 0.8 fm, as motivated in the 
previous sections. Secondly, together with the sources O", eq. (|2.5p . and C", which 
correspond to wo(x) = 1, we consider the following basis of wave- functions 

(4.7) 



Wi(x) = 


= iVri ^^,(|x-nL|), t = 

nGZ3 


: 1,2,3, 




wi(r) = 


= e"''/'^", ZJ2(r-) = re-''/'^% 


^■i{r) = 


^ g-r/(2ao) 



in order to construct the external operators. Here, we keep the physical length scale ao 
fixed in units of L by choosing oq = L/6 and the (dimensionful) coefficients Ni are set to 
normalize the wave function via a^ ^^ ^li'^) — ^^ ■ The same set of interpolating fields, 
for the same L, has been used in ref. [29] to determine the improvement coefficient ca 
by requiring the quark mass derived from the PCAC Ward identity to stay the same 
as the external states are changed. We therefore choose the already determined [29] 
optimal wave-function 

3 

u;^(o) = ^r/f ^cJi , r/(°) = (0.5172, 0.6023, 0.6081) . (4.8) 

i=l 

It suppresses the contribution of the ffist excited state in the pseudoscalar channel to 
the correlation functions under consideration (see section 2 in ref. [29]). 

The final result will turn out to differ significantly from the determination in ref. [10] 
at the two largest couplings only. We therefore did not recompute Za for small couplings 
and rather use the old estimates. 

4.1 Results 

We have two dynamical flavors of non-perturbatively improved Wilson quarks [30] and 
the plaquette gauge action. The improvement coefficients Cgw , ca were set to their non- 
perturbative values [29,30]. We chose T = 3/2L with periodic boundary conditions 
{6 = 0) in space and vanishing background field. We used the HMC algorithm with two 



pseudo-fermion fields as proposed in refs. [1,31]. The particular implementation has 
been discussed and tested in refs. [32] and [33]. Following the last reference we chose a 
trajectory length of r = 2 except for at /3 = 5.2 where we set r = 1. 
The normalization factor Za is given by [10] 



ZA{gl) 



lim v7iH 

m— »0 



/iA(2r/3, r/3, c^) - 2m/iA(2r/3, T/3, cu) 



-1/2 



(4.9) 



The definition of /p^(2T/3, T/3,u;) and the PCAC quark mass m can be found in 
ref . [10] . After performing the Wick contractions one realizes that disconnected quark di- 
agrams, where no propagator connects the two boundaries, contribute to /xy(3;o, 2/0, ^)- 
These can be shown to vanish in the continuum massless limit as a consequence of the 
conservation of the axial current [10]. In an improved theory they therefore amount 
to O(a^) effects on Zp^. By dropping them one obtains an alternative definition of Za, 
denoted Z'^^. With the kinematics of ref. [10] the difference between Z^ and Z'^'^ for 
/3 = 6/5q < 5.5 was found to be rather large, though consistent with O(a^) scaling. As it 
will become clear in the following this effect is very small in the computation presented 
here, which therefore significantly improves on the result in ref. [10]. 

In order to ensure a smooth dependence of Za on the bare coupling g^ and the 
correct scaling of discretization errors proportional to a^, we impose our normalization 
condition on a line of constant physics. This requires keeping all length scales fixed 
as Qq is varied. The lattice size L is set to approximately 1.8L* with L* given by the 
condition (f'{L*) = 5.5, where g'^{L) is the Schrodinger functional coupling defined in 
refs. [11,24]. The relation between — and go could be taken from ref. [34]. 

The bare parameters of our simulations and the results for Za are collected in 
Table [TJ Due to algorithmic instabilities caused by the appearance of very small. 











UJ = 


= Wo 


UJ = 


^7r(0) 


L/a 


/? 


K 


am 


Za 


^con 


Za 


rycon 
^A 


12 


5.2 


0.1355 


0.02121(36) 


0.788(15) 


0.784(4) 


0.784(16) 


0.7874(35) 


12 


5.2 


0.1357 


0.01434(48) 


0.775(10) 


0.769(4) 


0.777(11) 


0.7703(36) 


12 


5.2 


0.1358 


0.00907(39) 


0.776( 8) 


0.777(4) 


0.776(9) 


0.7788(37) 


12 


5.2 




^0 


0.766(18) 


0.773(15) 


0.769(20) 


0.774(16) 


16 


5.4 


0.136645 


0.00062(26) 






0.779(5) 


0.793(5) 


24 


5.7 


0.136704 


0.00072(14) 






0.808(5) 


0.802(3) 



Table 1: Bare parameters and simulation results (statistical errors only). At /3 = 5.2 the value of the 
renormalization constant extrapolated to the chiral limit is reported in the fourth line including the 
associated systematic uncertainty (see text for details). The number of decorrelated measurements used 
varies between 1200 (at (3 = 5.2) and 200 (at (3 = 5.7j. 

unphysical, eigenvalues in the spectrum of the Wilson-Dirac (SF) operator [3,35], at the 
coarsest lattice spacing we could simulate down to bare quark masses of about am k, 



0.795 

0.79 
0.785 

0.78 
^:^ 0.775 ^-4 

0.77 
0.765 

0.76 
0.755 



1 1 

only statistical error 
with systematic error 



0.005 



0.01 



0.015 



am 



0.02 



0.025 



Figure 3: Estimate of Z]^'^ at the chiral point for (5 = 5.2. The error-bar of the Riled circle at am = 
shows the statistical error whereas the error on the open circle includes our estimate of the 0{am) 
systematic uncertainty. 



0.01 only. Our estimates for the massless limits at (3 = 5.2 are just weighted averages of 
the numbers at the two lightest quark masses. The errors in Table [1] include a systematic 
uncertainty for possible 0{am) contaminations given by the difference between the 
determination at the heaviest mass and the described estimate. This uncertainty is 
added linearly to the statistical error. A similar 0{am) uncertainty would be obtained 
by comparing to a linear fit in am with all three points. Figure [3] illustrates the procedure 
for Z^°" and to = uj^(o)- It is clear from the plot that the dependence of Za on the 
quark mass is rather mild as expected for the "massive" definition of the renormalization 
constant [10] used here. Similar remarks obviously apply to the case uj = ujq, as Tabled] 
shows that the two wave- functions give consistent results for all values of n. At the other 
values of the bare coupling we simulated at very small quark masses. Given also the 
flat dependence observed at /3 = 5.2 we did not need to estimate an effect of order am. 
As anticipated, already at (3 = 5.7 the present result nicely agrees with the numbers in 
ref. [10]. Conversely, at /3 = 5.2 the new determination of Za is about 8% larger than 
the old one. These comparisons are summarized in Figure HI 

In view of the above discussion we smoothly interpolate the data for Z™"^ and 
CO = u;^(o)in the region 5.2 < /3 < 5.7 and those in ref. [10] for 7.2 < (3 < 9.6 (filled 
circles in Figured]) by the formula 



1 - 0.116ffg + 0.0115^ - 0.07256 , 



(4.10) 



where the coefficient of the term linear in g^ is fixed by 1-loop perturbation theory [36] 
and the last two coefficients are the result of a fit. We ascribe an absolute error to Za, 



N 



0.95 




0.75 - 



Figure 4: Comparison of the new results for ZX°" with the old ones for Za and Z^"^ in ref. [10]. The 
solid curve represents the interpolation formula in eq. i)4.iOI) . 



which decreases from 0.016 at /? = 5.2 over 0.005 at (3 = 5.4 to 0.003 for /? > 5.7. We 
note that our systematic 0{am) error at /3 = 5.2 is rather conservative, see Fig. El 



4.2 Other renormalization conditions of the Eixial current 

In order to get a further impression about residual cutoff effects, we also studied an 
alternative definition of Za- It is obtained in the same framework by replacing a light 
quark with a static one in the external operators O and C. By remaining with the flat 
wave-function and by denoting the static quark field on the boundary xq = by Ch we 
write 



""-h^^ 
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Of = ^ > : Ch(x)75C.(y) , or = j^Y. Ckx)75CUy) , (4.ii) 

x,y x,y 



where i = 1,2 is a flavour index. The flavour contractions in /^y are changed corre- 



hl,I 



spondingly. The correlator /^J^ for example is written 



M.lf N 

/AA(^o,yo) 






Y,e'''^{0't\AiTo{x) (A,t{y)\{r%,Of) (4.12) 



x,y 



and the massless normalization condition becomes 

72 A\,lt^rrin rrin\ _ i-hl 



zi/^;^(2r/3,r/3)=/r + o(a^ 



(4.13) 
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with 



n 



hi 



1 



/hl/nhl\ 



^ {oro, 



(4.14) 



As the fields X and Y in the correlator do not contain static fields, it is clear that dis- 
connected diagrams cannot appear and the static quark propagates from one boundary 
to the other. Static quarks are discretized through the HYPl and HYP2 static-quark 
actions [37], which have a relatively good signal to noise ratio in static-light correlation 
functions at large time separations. 



L/a 


T/a 


/5 


K 


am 


9 


Za [HYPl] 


Za [HYP2] 


12 


18 


5.2 


0.1355 


0.02121(36) 





0.7904(66) 


0.7882(55) 


12 


18 


5.2 


0.1357 


0.01434(48) 





0.7729(79) 


0.7734(73) 


12 


18 


5.2 


0.1358 


0.00907(39) 





0.7731(75) 


0.7736(70) 


8 


18 


5.2 


0.1357 


0.00662(75) 


0.5 


0.8284(80) 


0.8285(66) 


8 


18 


5.2 


0.1358 


0.00325(82) 


0.5 


0.8252(94) 


0.8206(82) 



Table 2: Simulations parameters and results for the alternative "static" definition of Za- 



The check was performed at the largest lattice spacing where ambiguities are ex- 
pected to be most pronounced. Simulation parameters and results are collected in 
Table [21 The results at L ~ 1.2 fm agree with those in the previous section, indicating 
again small overall cutoff effects. The values at L ~ 0.8 fm instead again suggest that 
the determination in ref. [10] suffers from large a^ effects, as both Za and Z™° there 
differ significantly (and especially for small quark masses) from the numbers in Table [21 



4.3 Renormalization of the vector current 

We also recomputed the renormalization constant Zy of the vector current in the new 
kinematics. Since it changes by less than 2% compared to ref. [10] at the largest lattice 
spacing, there is no reason to publish a new determination. 

5 Conclusions 

We have discussed possible kinematical enhancements of cutoff effects in lattice gauge 
theory determinations of QCD correlation functions. In this respect, the typical SF 
correlation functions improve significantly over those of composite local fields with pe- 
riodic boundary conditions (or large volume) . The origin of this difference is simply the 
difference in mass dimensions of the Schrodinger functional boundary fields compared 
to the usual local composite fields. Since the latter is at least three, time-slice corre- 
lators diverge at least as Xq at short distances. Such a steep behavior is difficult to 
approximate in a discretized theory. 

In the Schrodinger functional we identified a parametrically much weaker, but still 
relevant, kinematical enhancement of cutoff effects. It appears in the transition region 
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between approximately perturbative behaviour and dominantly non-perturbative one. 
Numerically we find that it appears e.g. for a T x L^ geometry with L ~ 0.8 fm , T « 
1.8 fm. In this kinematical situation a few hadronic intermediate states are relevant for 
the correlators and can produce a relatively steep decay of correlation functions even at 
a small quark mass. 

This means that the transition region from approximately perturbative to strongly 
non-perturbative is a relatively difficult one for numerical simulations. Discretization 
errors have to be investigated carefully. Fortunately we also saw that this region is 
rather narrow. With the step scaling method [19, 20, 22] it is typically bridged by one 
step. 

Avoiding the difficult region in the renormalization condition for the relativistic 
axial current, we have finally presented a significant improvement of the previous de- 
termination of Za [10]. The difference to the old one is of order a?, but it is up to 10 % 
at the largest lattice spacing considered. In our new determination we find that discon- 
nected contributions, which can be shown to vanish in the continuum limit, are very 
small - in contrast to the previous kinematical setup [10]. We also see a nice agreement 
with a renormalization condition where a static quark is present as a spectator in the 
Ward identity. The previous determination of Zy is confirmed. 
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